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1 Introduction 



The Nambu-Jona-Lasinio model [1] and generalizations thereof [2,3,4] have been used 
extensively to study the properties of mesons in free space and at finite temperatures 
and densities. It is an effective field theory of non-linearly interacting quarks which ex- 
hibits spontaneous and explicit dynamical chiral symmetry breaking. In the case of three 
flavours, it is mandatory to incorporate the 't Hooft six-fermion interaction to describe 
the breaking of the axial U(l) symmetry [3]. Mesons are bound quark- antiquark pairs 
in this approach and their properties can readily be calculated by solving the pertinent 
Bethe-Salpeter equations. However, no systematic study of three-point functions like the 
strong and radiative meson decays, has been performed so far. Previous attempts were 
limited to the leading term in the momentum expansion of the underlying quark-meson 
vertices [5] and thus do not account for the full dynamics of the model. Furthermore, 
these decays are also a good testing ground to find out the limitations of the model. This 
will be discussed in some detail later on. 

In what follows, we will work in flavour 577(3) and use the following Lagrangian 

C= G^X^f + (^ 75 A^) 2 ] 

+ K[det{^(l+7 5 )^} + det{^(l-7 5 )V'}] (1) 

where the flavour index % runs from to 8 with A = ^2/31, the A a are color matrices 
(a = 1, 8). As it stands, the Lagrangian is characterized by a few parameters: The two 
four-fermion coupling constants G 1 and G 2 , the six-fermion coupling K and the cutoff A, 
which is necessary to regularize the divergences. We will use a covariant four-momentum 
cutoff A = 1 GeV. Furthermore, to account for the explicit symmetry breaking, a quark 
mass term has to be added. We work in the isospin limit m u = m d and will use the 
current quark masses to fit the meson spectrum. Clearly, the coupling G\ is related to the 
properties of the pseudoscalar Goldstone bosons, G 2 can be fixed from the p-meson mass 
and K is necessary to give the rji]' mass splitting. This completely specifies the model 
and we are now at a point to consider its dynamical content. 

2 Formalism 

The basic object to consider is the triangle diagram which describes the coupling of the 
decaying meson (Mi) into the other mesons (M 2 ^ 3 ) or another meson (M 2 ) and a photon 7 
or two photons 71, 72 via the quark loop. Let us first consider the strong decays. Dropping 
all prefactors, the transition amplitude for the process Mi — ► M2M3 can be evaluated by 
working out (cf. fig.l) 

T(Mi -> M 2 M 3 ) = Tr(T Ml S F T M2 S F T M3 S F ) (2) 
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where T M . gives the i th meson-quark-antiquark vertex and Sp the propagator of the con- 
stituent quarks. The latter follows from minimizing the effective potential to one loop. 
The Bethe-Salpeter vertex functions relevant for our considerations are of scalar, pseu- 
doscalar and vector type 

T s = g s l®I 

T P = g P (l + hp j6) 75 <g> I (3) 
Ty = gvlv ® I 

Here, / is a generic symbol for the isospin structure and we have introduced scalar, 
pseudoscalar and vector Mgg-couplings. The coupling h p stems from the pseudoscalar- 
axial vector meson mixing. This is discussed in more detal in refs. [2,4,6,7]. Since we 
work in the isospin limit, no scalar-vector mixing arises in the SU(2) subgroup. In the 
scalar and pseudoscalar channels, a further complication is due to the AoA§ mixing which 
has already been discussed in ref.[8] in some detail. The solution of the corresponding 
Bethe-Salpeter equations is standard and we do not exhibit any details here. 

Let us briefly elaborate on the connection between the various transition amplitudes 
and meson-meson coupling constants. Consider first the decay of a scalar (S) into two 
pseudoscalars (P). The transition amplitude is a purely scalar function, called Ts, and 
we have 

\Pc\\Ts\ 2 



T(S -> PP) 



8ttEI 



GIpp = ^ (4) 

with \p c \ the momentum of an outgoing particle in the rest frame of the decaying particle, 
\p c \ = ((s — (mi + m 2 ) 2 )(s — (mi — m2) 2 )/4s) 1 / 2 and s = m 2 (Mi). For the decay of a 
vector (V) into two pseudoscalars, one has 

T^V^PP) = ( Pl -p 2 ),G VP p 
It? \ 3 n 2 

r(v ^ PP) = m u \ pp (5) 

Qirniy 

Finally, for the reaction V — > VP we find 

T^V ^ VP) = e^pfflP 

T(V ^VP)J^f^- (6) 

In the case of the radiative decays, one can use the same formalism since the photon 
behaves much like a vector particle. Since our Lagrangian contains no tensor interaction 
term, the pertinent photon-quark-antiquark vertex takes the minimal form 

r 7 = |(A 3 + ^A 8 ) (7) 

For the decay of a pseudoscalar into two photons, we have a structure similar to the one 
for V — > VP, the only difference being that in the formula for the width T(P — > 77) one 
has a factor 1/2 instead of 1/3 which is the reduction in plolarization degrees of freedom 
for a massless particle. 
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3 Estimates of higher order effects 



The description of mesons as qq pairs has been proven to be quite successful in the 
calculation of the meson mass spectrum. However it may not account for other properties 
of some mesonic resonances, such as their decays. We anticipate that this is indeed the 
case for the strong decays p — > nn and K* — > Kir, which come out very small in the 
one loop calculation. The insufficiency of a qq description of the mesons was already 
emphasized by Krewald et al [9] in the context of the pion electromagnetic form factor. 

Due to the failure of the one loop approximation for the calculation of those decays, 
it is necessary to estimate the magnitude of higher order effects, such as the two loop 
corrections. The full calculation is, for the moment, out of the scope of the present work. 
A simple estimate of such effects for the mesonic decays can be obtained by calculating a 
dressed meson propagator, as shown in fig. 2. The dashed loop refers to nn or Kn states, 
in the case of the p and K* propagators, respectively. By bare propagator we denote the 
meson described as a qq state. 

The full propagator reads 

G a fi = G° a(3 + G° aX S Am G^p (8) 

where the bare propagator is 

n Q _ 9af3 



and Sa m is the meson loop given by 



— m 



v 



q 2 — m v 



(9) 



S AM = 8G 2 VPP f ^-((2k-q) x )((2k-q),)S(k,m P1 )S(k-q,m P2 ) (10) 
JAi {2n)^ 

with S(k,mpi) the propagators of the mesons obtained in the one loop order, m P \ and m P2 
being their masses. The remaining factors in the integrand correspond to the coupling of 
the vector mesons to the pseudoscalars in the first order calculation (5). The covariant 
cutoff A 1 needs not to be the same as A. One obtains: 

Ga " = q 2 -m\-C 

with C = C ( G VPP ,q 2 ,mp\ ,mp 2 , Ai ). The quantity C has a cut for q 2 > {m P i + 
m P2 ) 2 and one rewrites G a p as 

G af 3 = T af3 — rn r-f — ~ 7 z 2 w (12) 

q 2 — m v — i Im C (m v ) Z 

provided that the renormalization factor Z is roughly constant around the physical mass 
my and where the renormalization factor is: 



d_ 

dq 

The decay width of the vector meson is finally 



q 2 =m v 



(13) 



4 



r ( v pp) = hC W» z . (M) 

ffly 

Using this scheme the amplitudes for radiative decays of p and if* have then to be also 
multiplied by \fZ (wave function renormalization). 



4 Results and discussion 

At the one loop level we use the meson spectrum and decay constants to fix the param- 
eters of the model. For A = 1 GeV, GiA 2 = 3.95, G 2 A 2 = 5.43, KA 5 = 42, m u = m d = 4 
MeV and m s = 115 MeV, we find the following meson masses (the experimental values 
are given in parentheses for comparison): M n = 136.5(139.6), Mk = 497.5(497.7), M v = 
549 (548.8), M v , = 936 (957.5), M p = 775 (768.3), M u = 764 (781.95), M K * = 898 (891.6), 
M<j> = 990(1019.41) and M ao = 970(983.3) (all in MeV). For the pion and kaon decay 
constants, we have F w = 93.9 MeV and F K = 96.6 MeV, i.e. the ratio F K /F W is too small, 
which is a common feature in this kind of models. We find an overall satisfactory descrip- 
tion of the meson spectrum together with reasonable values for the vacuum expectation 
values of the scalar quark densities uu and ss, - < uu > 1 / 3 = 272 MeV (225 ± 35) and 
< ss > I < uu >= 0.74(0.8 ± 0.2). 

In table 1, we show the results for the strong decays in comparison to the empirical 
values. Obviously, for states in the quark-antiquark continuum the results are not reliable 
as indicated by the decay $ — > up. Also, for our set of parameters the decay $ — > KK is 
kinematically forbidden. The large width of the $ — > np decay is due to the too strong 
flavour mixing induced by the six-fermion interaction proportional to K. This could 
presumably be cured by including more terms in the Lagrangian like e.g. in ref. [6]. 

As for the strong decays of p and K* the one loop order calculation is clearly insufficient 
to account for the respective empirical widths. Using the simple approximation scheme, 
described in the previous section, to include the second order effects and approximating 
the meson propagators in the loop by propagators of structureless particles, the results 
improve by about a factor of 2. We think, therefore, that it is mandatory to consider a 
more complex multiquark structure for the p and K* mesons. The strong p decay width 
including second order effects is still quite small as compared to the experimental one, but 
this number should be understood only as a guide for the order of magnitude of higher 
loop corrections. We notice that the parameters could have been adjusted in order to 
have the correct decay width for the p and the KSFR relation fulfilled, but at the cost of 
having bad values for the radiative decays. 

Let us now turn to the radiative decays. In table 2, our results are summarized. We 
find an overall satisfactory description of the data, the main exceptions being the widths 
for $ — » 777 and for the K* — > Kit. In the first case this is, again, an artifact of the 
interaction Lagrangian used and also, since the $ lies in the unphysical quark-antiquark 
continuum, should not be considered significantly troublesome. The fact that the ratios 
Typp/Typp and ^vp-y/^vp-y are l &r g er f° r the K* decays than for the p decays (both 
strong and radiative) is consistent with the small ratio Fk/F^. We have also investigated 
the case K = (no six-fermion term). Reducing the strange quark mass to 80 MeV [2], 
which is necessary to find a decent fit to the spectrum, one is not able to get a satisfactory 
description of strong and radiative decay widths. After finishing this work, we became 
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aware of a preprint by Takizawa and Krewald [10], who deal with the radiative decays 
7t° — > 77 and rj — * 77 in a similar model. Their Lagrangian contains the four-fermion 
interaction proportional to G\ and the six-fermion determinantal interaction. While 
their results are similar to ours, we disagree with their conclusions at various places. 
First, in the case of tt° — > 77 they remove the cut-off to find agreement with the current 
algebra prediction. This is, however, not a consistent procedure since once the cut-off 
is introduced, the effective theory is defined and should not be altered in the process of 
calculating various quantities. Second, for calculating the width of 77 — > 77, they use the 
empirical rj mass, which is 17 per cent larger than the value they find within the model. 
This, of course, alters substantially the result for this particular width. Comparing their 
results with ours, we also find a satisfactory description for these two particular radiative 
decays. It should be obvious, however, that the model is somewhat too crude too draw as 
far reaching conclusions as done in ref. [10]. As long as one is not able to properly account 
for the SU(3) breaking effects in the pseudoscalar decay constants, it is doubtful that one 
can make firm quantitative statements about such breaking effects in other processes. 

In summary, we have used the generalized three-flavour NJL model to calculate strong 
and radiative meson decays (three-point functions) taking the full solution to the Bethe- 
Salpeter equations in the one-loop approximation. A reasonable agreement to the ex- 
perimental data is obtained, with the exception of the p and K* decays. This seems to 
indicate that a more complex multiquark structure should be accounted for these mesons. 
This conjecture is supported by a simple estimate of the two loop order corrections. We 
also demonstrated the importance of the flavour- mixing determinantal six-fermion in- 
teraction. Further studies including also the effects of isospin-breaking seem necessary 
to understand some fine details of the mesonic interactions at low energies and to gain 
insight into effects of SU (3) breaking on various observables. 
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TABLES 





p — > TTTT 


K*+ -> 7i+K° 


K*+ -> 7T°^+ 


a -> tt?7 


— > 7rp 


NJL(I) 


52.0 


18.0 


9.0 


74.5 


1.5 


NJL(II) 


94.0 


38.2 


19.1 






Exp. 


151.5 ± 1.2 


38.6 ±0.6 


16.7 ±0.3 


57± 11 


0.6 ±0.3 



Table 1: Strong meson widths for various decays in units of MeV: (I) calculated in one 
loop order; (II) with estimates of two loop order included. 





7T° — > 77 


77 -»• 77 


p 1 * 1 — > 7T ± 7 


P -»• ^7 


NJL 


7.9 • 10" 3 


0.77 


60.1 


60.4 


Exp. 


(7.7 ±0.6) • 10" 3 


0.46 ±0.04 


67.1 ±7.6 


57.6 ± 10.7 




UJ — > 7T°7 


c<j — > r?7 


K* + -> K+7 


$ — > r?7 


NJL 


762 


6.3 


92.0 


259.0 


Exp. 


716.6 ±43.0 


4.0 ± 1.9 


50.3 ±4.6 


56.7 ±2.8 



Table 2: Anomalous and radiative meson decay widths in units of keV, calculated in the 
one loop order. The second order corrected decays are r p7r7 = 63 keV and T K * Kl = 98.9 
keV (Ax = 1.3 GeV). 



FIGURE CAPTIONS 

Fig. 1: Quark triangle diagram to calculate the strong meson decays. In the case 
of radiative decays, one has to substitute the third meson M 3 by a photon, and for the 
anomalous decays M 2 and M 3 by two photons. 

Fig. 2: The dressed vector meson propagator (thick line) includes a 2-pseudoscalar 
excitation (dashed line), which is a 2nd order effect. The one loop order (qq state) is 
denoted by the thin line, the bare meson propagator. 
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